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We determined the eigenstates of a single electron in a 
parabolic anisotropic 2D quantum dot in a magnetic field. 
Using obtained expressions for these eigenstates, we study 
the spin coupling J between two electrons located in two lat- 
erally coupled anisotropic quantum dots (QD). The exhange 
coupling J is calculated using the Heitler-London and Hund- 
MuUiken approaches. We found that the exhange J changes 
sign at certain values of parameters of the system, in partic- 
ular at certain anisotropy of the QD. Therefore, we present a 
new method to switch on and off the spin coupling between 
QD: switching by means of changing shape of the QD. Switch- 
ing the spin coupling is essential for quantum computation 
using electronic spins as qubits. We note that our calcula- 
tions can be applied to the system of vertically coupled QD 
as well. 



I. INTRODUCTION. 

The last decade has seen a great interest in quantum 
dots, i.e. little islands of electron gas confined by an 
artificial potential (see Ref. One of the most chal- 

lenging proposals concerning QD is the idea of using the 
electron spin in quantum dots as the basic information 
carrier (the qubit) in a quantum computer (see, for exam- 
ple, Ref. Ql). Quantum logic gates between these qubits 
are effected by allowing the electrons to tunnel between 
two coupled quantum dots. For the application of cou- 
pled quantum dots as a quantum gate, it is important 
that the coupling between the spins can be switched on 
and off via externally controlled parameters such as mag- 
netic field and interdot distance. In the recent paper of 
G.Burcard et al.^ the spin interaction for two laterally 
coupled isotropic QD defined in a two-dimensional elec- 
tron system (2DES) was calculated as a function of these 
external parameters and it was shown that the interac- 
tion J can be switched on and off with exponential sen- 
sivity to that parameters. In the present paper, we con- 
sider a difi^erent setup consisting of two laterally coupled 
anisotropic quantum dots with a magnetic field applied 
perpendicular to that 2DES (see Fig.l). We especially 
interested in calculation of interaction J as a function of 
anisotropy of the QD since our proposal is to switch on 
and off the interaction J by changing shape of the QD, 
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e.g. by applying gate voltage to each quantum dot. We 
note that in the setup being proposed in Ref. j^] it was as- 
sumed that the interdot distance is controlled by varying 
the barrier height between the QD. In real experimental 
setup varying the barrier height might lead to changing 
shape of the QD from circular to elliptic. Also, probably 
in future quantum computers based on large integrated 
circuts of QD it will be simplier to change shape of the 
QD (and, in case of neceserrity, interdot distance), than 
change interdot distance alone, preserving circular shape 
of the QD. This is one of the reasons which motivated 
our study. 

In order to study coupled anisotropic quantum dots, 
one needs to find firstly the eigenstates of a single elec- 
tron in a single anisotropic quantum dot in a magnetic 
field. The latter problem is of fundamental interest in 
mesoscopic physics. Surprisingly, it seems that this im- 
portant problem was not solved so far. 

The eigenstates of a single two-dimensional electron 
confined by a radial potential of the form imeW^r^ with 
a magnetic field perpendicular to the plane of the system 
are called Fock-Darwin states. That system was investi- 
gated by Fock'^ and Darwin^ a long time ago. The system 
where an electron is confined by an anisotropic potential 
of the form ^me(uj'^x'^ + cUyi/^) was studied only recently 
by A.V. Madhav and Tapash Chakraborty.^ Their stud- 
ies were motivated, in particular, by the experiments on 
anisotropic quantum dots (see Ref. Q). They obtained 
energy eigenvalues of the system, but not the eigenfunc- 
tions of it. In their paper, the transformation of coordi- 
nates and momenta of the system was performed. In the 
new Hamiltonian, variables are separated and one can 
easily obtain both the energy eigenvalues and the eigen- 
states of the system. However, in order to obtain the 
eigenstates of the system in old (initial) coordinates, one 
needs to transform the eigenfunctions of the new Hamil- 
tonian. It was not done in Ref. |^ . In the section II of the 
present paper, we performed such a transformation and 
obtained exact analytical expressions for the eigenstates 
of that system. Note that the system of an anisotropic 
quantum dot in a magnetic field was considered also in 
Ref. , where only approximate solution for the ground 
state of the system was obtained (see Appendix D in that 
paper). In the section HI we study coupled anisotropic 
QD using expressions obtained in section II. Following 
Ref. we employ Heitler-London approximation and 
Hund-MuUiken approach. In section IV we give brief 
conclusion. 
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II. SINGLE ANISOTROPIC QUANTUM DOT. 

Let us consider a single electron in a lateral anisotropic 
parabolic confinement potential in the presence of a per- 
pendicular magnetic field. The Hamiltonian is 
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(1) 



where rrie is the effective mass of the electron and the 
confinement potential is 



yconf{x,y) = Kn^{wlx^ +ujly'^). 



(2) 



Choosing the symmetric gauge vector potential A = 
lB{-y,x,0), we get 

H = -^[pI + ^\x^ +pI + + meUJciyPx - xpy)], 
Imp, " 

(3) 
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Following Ref. Q, we make the following transforma- 
tions: 



x = ^qi+ r]p2, 
y = m2 + VPi, 
Px = MPi - C92, 
Py MP2 - C9i, 



(4) 



where 
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The transformations are consistent with the com- 
mutation relations [pi,qj\ = —ifi5ij,[qi,qj\ = and the 
resulting Hamiltonian is diagonal (see also Ref. |^) : 



H = T^i^lpl + alpl + Plql + PhD 

where 



(6) 
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One can easily obtain eigenvalues and eigenstates of the 
Hamiltonian: 
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'0m,n(gi, 92) = Cm,„e ^ e 2 H™ (71(71 )H„ (72^2 ) , (9) 
where Cm,n = ( ^2^^+^Lin! )^^^' = a,(3^/me, jf = 

f3i/{aih), Hk is the Hermit polynomial. The question 
is how one can return back to coordinates (x,y), i.e. how 
to obtain eigenstates '^m,nix,y) of the Hamiltonian (||). 

Since ip{qi,q2) = (91,92!*), *(a;,?/) = (a;,y|*), and 
I*) ^ I-!^ d,qidq2\qi,q2)^j{qi,q2), we have the following 
formula for the eigenstates '^rn,n{x,y)'. 



^m,n{x,y) 



+ 00 



dqidq2 {x,y\qi, q2)^m,n {qi , 92) (10) 



Let us denote $(a:, y, gi, 92) — {x,y\qi,q2)- Function 
^{x,y,qi,q2) is the eigenstate of operators (71,92 in the 
(x,y) representation and should obey the following rela- 
tions: 



92$ = [^J,y - Tipx]'^ = 92$. 



(11) 



On the other hand, <i>*(x, y, 91, 92) = {qi,q2\x,y) ('*' de- 
notes complex conjugating) and therefore it should obey 
the following relations: 



= [fiqi + ?7P2]$* = a;$*, 
^ [M92 +7?P2]$* = y$*. 



(12) 



Therefore, we finally have the following system of dif- 
ferential equations defining <i>(a;, y, 91, 92): 

{ihri-§- + nx)<^{x,y,qi,q2) = 9i$(a;, y, 91, 92), 



{ihrig^ + ^iy)<^{x,y,ql,q2) = 92$(a;, y, 91, 92), 
{^i^'n^+ mi)^*{x,y,qi,q2) = x$*(x,?/,9i,92), 
m2)^*{x,y,qi,q2) = ?/$*(x,?/,9i,92). 



One can find that solution of this system has the fol- 
lowing form: 



"J'(a;,y,9i,92) = Ccxp | — (/^a;?/ + ^9i92 - qiy - 92a;) 

(13) 

where C is a constant. From ( p^p^p |) we have 

^m,n{x, y) = Cm^ne'^ Im.nillVX, J2'^y), (14) 

where Cm,n = C'c™^„/7i72, v = {fi'^ + ri'^h'^-fljD~'^/'^,and 

Im,n{zi,Z2) = 
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where a = /// (77 717172), d = Vl + cr^. Let us firstly find 
the ground state of our system \E'oo(x, y). DifRrentiating 
expression for 100(2:1,22) over zi and Z2 and solving ob- 
tained difhrcntial equations, one can find that 



Im{zi,Z2) = Cooe 



2 ^2 
f -^„-»0-ZiZ2 



(16) 



where Coo is a constant. Therefore, the ground state of 
the system is 
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where normalization constant is Coo — \J "^"^^ ; = f^7i j 
Kj, = 1^72. One can see that CCqo = Hxi^y/v. Let 
us now note that H„(g) = (2g)" - ^^^^^(2^)"^^ ^_ 

ri(»-l)(K-2)(rt-3) _ 

Let Lfc(2i,22) = /+^(2g2)'=i^(gi,g2,22)e-*«^^i'^dqidq2 
, where F is some function. One can see that 

2i dLk-i 
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A similar procedure for the coordinate Z2 can be exe- 
cuted. Now one can easily find the following expression 
for /m,„(zi, 22) : 



Im,n{zi, 22) = -D™I?"^/oo(2l, 22): 

where operators are defined as 

i d 
d da 

2iY d'' fc(fc-l) (2i\''-^ 
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So, we finally have the following expression for the eigen- 
statcs of the Hamiltonian (S): 
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where I(x,y) = exp{— -^^^^^ ^ iaKxHyXy}. The 

lowest excited states of the system are 
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(o-Kj^y - iKj;a::)i^oo(2^,2;), 

(an^x - iKyy)FQo{x,y), (22) 

[K2;Kya;j/(l - fT^) + 



*oi(2;,y) 
*io(a;,?/) 
*ii(a;,2/) 

where 

Foo(x,y) = e'^^i-'-'^e-^-^ = CoV^'oo(a;, y). 
Note that in degenerate case of isotropic dot we have 
rii = fl2, V = a = 1, and the eigenstates ( p2| , pT[ ) coincide 
with the corresponding Fock-Darwin states. 
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III. COUPLED ANISOTROPIC QUANTUM 
DOTS. 

A. ModeL 

The Hamiltonian which we use for the description of 
two laterally coupled quantum dots is 



h{r,p) 



1=1,2 
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p-£A(r))\nr), 



(23) 



C = 



K ri - r2 



where C is the Coulomb interaction and h is the single- 
particle Hamiltonian. The dielectric constant k and the 
effective mass TOg are material parameters. For the lat- 
eral confinement V we choose the quartic potential 



V{x,y)^ 



2\2 
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(24) 



which separates (for x around ±a) into two anisotropic 
harmonic wells in the limit of large inter-dot distance, i.e. 
for 2a 3> 2aBx,y , where a is half the distance between 
the centers of the dots, and the lateral effective Bohr radii 
o-Bx,y = \/'hlmeUix,y are a measure for the lateral exten- 
sion of the electron wave function in the dots (see Fig. 1). 
In Ref. 1^ the potential similar to (j2^) was considered, 
but there uj.j. = Wy = tuo- It was also shown in Ref. ||] 
that spin-orbit contribution and Zeeman splitting can be 
neglected in relevant cases. For material parameters, we 
choose that of GaAs {k = 13.1, nie = 0.067m, where m 
is the mass of the electron). 



B. Heitler-London approach. 

Here we consider firstly the Heitler-London approxi- 
mation, and than refine this approach by including (see 
subsection C) double occupancy in a Hund-MuUiken ap- 
proach. In the Heitler-London approach, one starts 
from single-dot ground-state coordinate wavefunctions 
</j(r) and combines them into the (anti-) symmetric two- 
particle state vector 



|12)±|21) 



(25) 



the positive (negative) sign corresponding to the spin sin- 
glet (triplet) state, and S = J d'^r(fX^{r)^_^{r) = (2|1) 
denoting the overlap of the right and left orbitals. A 
non-vanishing overlap implies that the electrons tunnel 
between the dots (see also Ref. |2j). Here, </J-Q(r) = (r|l) 
and (p+a{^) = (r|2) denote the one-particle orbitals cen- 
tered at r = (=Fa, 0), and — \i)\j) are two-particle 
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product states. The exchange energy is then obtained as 
J = et - = (*_|ir|*_) - {^+\H\'^+). The single-dot 
orbital for anisotropic harmonic confinement in two di- 
mensions in a perpendicular magnetic field is the ground 
state being obtained in section II. The ground state cen- 
tered at the origin is 



(26) 



where 7 — — t^^)/?7?i, and Ky were introduced in 
section II. 

Shifting the single particle orbitals to (±a, 0) in the 
presence of a magnetic field we obtain ^p±a{x,y) = 
exp(±i?/a/2/^)(^(a;=Fa, y). The phase factor involving the 
magnetic length Ib ~ yJfic/eB is due to the gauge trans- 
formation A±a = B{-y,x^a,Qi)/2 ^ A = B{-y, x,0)/2 
(see also Ref. Q). We rewrite the Hamiltonian, adding 
and subtracting the potential of the single left (right) dot 
for electron 1(2) in H. The Hamiltonian then takes the 
form H = hP_^{ri) + h°^_^{r2) + W + C, where /i±a(r.,) = 
{pi-eA{ri)/c)^ /2me+meUjl{xiTaf /2+meUjlyf /2, and 
W{ri,r2) = Ei=i,2 yixt,yt) - rueujliixi + af + ix2 - 
a)^)/2 — meLUy{y1 + y|)/2. The formal expression for J 
is now 



J = 



2^2 



;i2|C + T4^|12) 



Rc(12|C + iy|21) 
52 



(27) 



where the overlap becomes S = exp(— K^a — ). 
Here we denoted Lf = (07 — a/2lg)~'^. Evaluation of the 
matrix elements of C and W provides us with the result, 



J = 



(1 - kIk-^) loiKW" 



CK,eS''-"'K (1 - >iIk-^) Io(k^'/;') 



y 

,2 , -2\ 



/sinh(2«;2a2 ^.2K-2r2^^ ^28) 



where c = \f2pne^ j K is the complete elliptic integral 
of the first kind, and Iq is the zeroth order Bessel func- 
tion. The first and second term in Eq.(28) are due to 
Coulomb interaction C, where the exchange term enters 
with a minus sign. The last term comes from the confine- 
ment potential W . Note that in the degenerate case of 
isotropic QD our result coincide with that being obtained 
in Ref. §. 

We plotted J as a function of lOy and lOx at different 
values of magnetic field B in Fig. 2. Note singlet-triplet 
crossing (J = 0) at certain parameters of the system. 



C. Hund-Mulliken approach. 

Following Ref. |^], we now apply the Hund-Mulliken 
approach to calculate the exchange energy of the 



double-dot system. We introduce the orthonormal- 
ized one-particle wave functions $±a = {}P± ~ 
.9¥'Ta)/\/l ~ 25(7 + where S again denotes the over- 
lap of (p^a with (y9+a and g — (1 — \/l — S'^)/S. Us- 
ing $±aj we generate four basis functions to which we 
diagonalize the two-particle Hamiltonian H: the states 
with double occupation, 5'^^(ri,r2) = <i>±a(ri)'i'±a(i"2) 
and the states with single occupation, \E'j_(ri,r2) = 

[$+a(ri)$_a(r2) ± $_a(ri)$+a(r2)]/V2. Diagonaliza- 
tion of 



H = 2e 
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2e + L/h/2 
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V+ ± 



yields the eigenvalues es± 
^/U'^/A + Ai^, €so ^ 2e + UH~2X + V+ (singlet), and 
Et = 2e + V- (triplet), where the following quantities were 
introduced : 

V = V--V+ = (*L|C|*1) - (^''+|C|1'"^), (30) 

x = {^i^c\^%:), 

Uu = U-V++X 

The exchange energy is the gap between the lowest singlet 
and the triplet state 

j^,,-e,_=V~^ + yUl + ml. (31) 

The results of evaluation of matrix elements are given in 
Appendix. We plotted resulted J as a function of ojy and 
LUx at different values of magnetic field B in Fig. 3. 



IV. CONCLUDING REMARKS. 

In summary, we obtained the eigenstates of two- 
dimensional anisotropic oscillator in a magnetic field (Eq. 
(pT|)) which in degenerate isotropic case become the well- 
known Fock-Darwin states. We used obtained analytical 
expressions to study a system of two coupled elliptic QD 
in a magnetic field. We have calculated the exchange 
energy J{B,a,uJx,ujy) between spins of coupled QD as 
a function of magnetic field, interdot distance and pa- 
rameters uJx,ijJy defining shape of a single quantum dot. 
We have shown that by varying Ux or Uy the exchange 
coupling J can be switched on and off. This opens up 
the alternative way (compared to that being investigated 
in Ref. ) of performing quantum computing operations 
using coupled QD as quantum gates. In the present pa- 
per, we didn't consider spin-orbit contribution, Zeeman 
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splitting and interplay between nonadiabaticity of vary- 
ing the system parameters and excitation of higher en- 
ergy levels. We consider these questions in a separate 
paper (see Ref. (H). We emphasize that our results can 
easily be applied to the system of two vertically coupled 
three-dimensional QD with in-plane magnetic field B be- 
ing studied in Ref. |^ , since the Hamiltonian of that prob- 
lem separates on _B-independent part, which is merely a 
harmonic oscillator, and i?-dependent part which actu- 
ally coincides with the Hamiltonian of two laterally cou- 
pled anisotropic QD being studied in the present paper. 

To conclude, we believe that our results bring closer 
the day the first quantum computer^~^° will begin to 
work. 



Fi = K^K (1 - nlnf) , (A8) 
F2 = KyK (1 - kIk~^) e'^'^'h (tila^) , (A9) 
F3 = (1 - KlKy^) eN ''"lo '^;') , (AlO) 

F4 = Kj;K (1 - kIk'^) e " "'4 ~ X 

(All) 

where I„ denotes the Bessel function of n-th order. 
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APPENDIX A: HUND-MULLIKEN MATRIX 
ELEMENTS 

Here, we give the explicit expressions for the matrix el- 
ements defined in Eqs. ( ^9| ) and (^0|). The single-particle 
matrix elements are given by 



3 m.ujl 3 52 



32 nta^ 



81-52 kI 



(1 + nla^) 



2mp 



2me 
3 S 



1 



where S = exp(— K^a^ 

The (two-particle) Coulomb matrix elements are for- 
mally equal to that given in Ref. ||^ , where only Fi have 
to be changed. We give here the complete set of expres- 
sions for convenience: 



52 
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^ (1 + Kla') 



(Al) 
(A2) 
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V+ = cN^ (4g2(l + ^2)^^ ^ ^ ^2^2p^ 

+ 4g^F3 - leg^Fi) , 
V-=cN\l-g'f{F2-S'F3), 
U = cN'^ ((1 + / + 2g^S^)Fi + 2g^F2 

+ 2g^S^F3 - Sg^Fi) , 
X^cN^[{{l + g^)S^ 



2g^) Fi + 2g^F2 



2g^S^F3 - Sg^FJ 



(A3) 
(A4) 

(A5) 
(A6) 



w = cN^ {~g{l + g^)il + S^)Fi - g{l + g^)F2 
~ g{l + g')S^F3 + (1 + 6.92 + /)5F4) , (A7) 



with N = 1/Vl - 2Sg + g^, c 



f , and 5 = (1 
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Figure captions. 



Figure 1. Two coupled anisotropic QD with one valence 
electron per dot. Each electron is confined to the xy 
plane. The magnetic field B is perpendicular to the 
plane, i.e. along the z axis. The quartic potential V{x, y) 
is given in Eq. (|2^ ) and is used to model the coupling of 
two anisotropic harmonic wells centered at (±a,0, 0). 

Figure 2. Exchange energy J in units of meV plot- 
ted against anisotropy ujy/ujx, as obtained from Heitler- 
London approximation. Upper graph: Ux was fixed 
(hux — 3 meV), tUy was varied; interdot distance a— 0.7 
ttBx- Bottom graph: ujy was fixed {hojy — 3 meV), Ux 
was varied; interdot distance a is the same as in the up- 
per graph. For both graphs, solid line: magnetic field 
B=0.8T; dashed-dot line: B=1.2T; dotted line: B^iT. 

Figure 3. Exchange energy J in units of meV plot- 
ted against anisotropy coy/ujx, as obtained from Hund- 
MuUiken approach. Upper graph: cOx was fixed {hujx — 3 
meV), LUy was varied; interdot distance a— 0.7 asx. Bot- 
tom graph: LUy was fixed {hujy — 3 meV), ujx was varied; 
interdot distance a is the same as in the upper graph. 
Solid line: 5=0.8 T; dashed-dot line: 5=1.2 T; dotted 
line: B=4 T. 
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This figure "figl.gif" is available in "gif" format from: 
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